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A recent experimental claim of the detection of analogue Hawking radiation in an optical system 
[PRL 105 (2010) 203901] has led to some controversy [PRL 107 (2011) 149401, 149402]. While this 
experiment strongly suggests some form of particle creation from the quantum vacuum (and hence 
it is per se very interesting), it is also true that it seems difficult to completely explain all features 
of the observations by adopting the perspective of a Hawking-like mechanism for the radiation. For 
instance, the observed photons are emitted parallel to the optical horizon, and the relevant optical 
horizon is itself defined in an unusual manner by combining group and phase velocities. This raises 
the question: Is this really Hawking radiation, or some other form of quantum vacuum radiation? 
Naive estimates of the amount of quantum vacuum radiation generated due to the rapidly changing 
refractive index — sometimes called the dynamical Casimir effect — are not encouraging. However 
we feel that naive estimates could be misleading depending on the quantitative magnitude of two 
specific physical effects: "pulse steepening" and "pulse cresting" . Plausible bounds on the maximum 
size of these two effects results in estimates much closer to the experimental observations, and we 
argue that the dynamical Casimir effect is now worth additional investigation. 

PACS numbers: 42.65.Re, 04.70.Dy, 42.65.Hw 



I. BACKGROUND 

It has been theoretically demonstrated that electro- 
magnetic field modes propagating in a dielectric medium 
with a spacetime dependent refractive index can closely 
resemble field modes in a black hole background geome- 
try pQ . The idea was further developed (and analyzed in 
more realistic contexts) in the sequence of papers [2HZ]- 

The central idea of that work (and follow-up develop- 
ments) is that a sufficiently intense, localized and mov- 
ing refractive index perturbation (RIP) can give rise to 
a region in the interior of the perturbation where elec- 
tromagnetic modes cannot propagate in the direction of 
motion of the perturbation. 

That is, a mode slows down while climbing the trailing 
edge of the RIP, and if the RIP is sufficiently intense and 
fast, it is possible for the mode's group velocity to vanish 
in the frame of the RIP. 

The condition for the existence of such a "blocking 
region", associated with a RIP of peak refractive index 
n mnx = n o + V traveling at velocity vpjp on top of a 
background refractive index uq, would be 
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This can equivalently be rephrased as a rather tight con- 
straint on the refractive index 
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Recently, in Ref. [8j Belgiorno et al. reported an inge- 
nious experimental realization inspired by this proposal, 
and it was extensively argued that the observed photon 



emission could be interpreted as the analogue of spon- 
taneous quantum Hawking radiation. Additional exper- 
imental and theoretical details are available in Ref. [9], 
and some critical analysis and response is available in 
Refs. PHI E] • Specifically, in this experiment an ultra- 
short laser pulse was used to create a traveling RIP in a 
transparent dielectric medium (fused silica glass) through 
the so-called non-linear Kerr effect (see e.g. [T]). The au- 
thors of Refs. [HIS] report experimental evidence of unpo- 
larized photon emission, for which they suggest the most 
compelling explanation is a Hawking-like radiation due 
to the presence of the above mentioned blocking region 
which simulates a black hole-white hole system. 

In closing this section, let us stress that this experi- 
ment should be contextualized in the wider framework 
of analogue models of gravity [H] and that several other 
experiments have been recently performed (or planned) 
in this context. To cite a few, the classical stimulated 
Hawking emission experiment of Weinfurtner et al. [13] , 
the hydraulic jump white hole of Jannes et al. |14) . and 
the negative phase velocity experiment of Rousseaux et 
al. [15]. There is also a significant quantity of theoreti- 
cal work on these more general analogue spacetimes, and 
their implications regarding the Hawking effect [TBHH]- 

Overall, we consider this RIP-based experimental sce- 
nario to be a very promising avenue in the ongoing search 
for a direct observation of some analogue of quantum vac- 
uum radiation in condensed matter systems, and we hope 
that the analysis here will shed some light in this direc- 
tion. 
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II. EXPERIMENTAL SITUATION 

The most salient features of the experimental observa- 
tion are the following: The fused silica was illuminated 
with a laser pulse of 1 ps timescale, and the input energy 
was varied in the 100-1200 [iJ range. Radiation from the 
filament was then collected at 90 degrees with respect to 
the laser pulse propagation axis. This arrangement was 
chosen in order to strongly suppress, or eliminate, known 
spurious effects (e.g. Cherenkov radiation, fluorescence, 
etc.) [8]. The observed emission spectra were reason- 
ably well fitted by Gaussian profiles. These Gaussian 
fits showed an increasing peak wavelength of emission, 
and increasing overall flux, with increasing pulse inten- 
sity. Characteristic numbers for the observations were 
the following: the emission is centered around a peak 
wavelength of about (850 ± 25) nm (see Fig. 3 of [8,) 
and has a bandwidth that appears to depend on the pulse 
intensity. While the determination of the actual size of 
such a bandwidth is fitting dependent, it is evident that 
the most of the photon count is concentrated in the range 
from 700 to 1000 nm. 

The interpretation of the observed spectrum as an 
optical analogue of spontaneous quantum Hawking ra- 
diation is based mainly on the above mentioned win- 
dow condition. In fact, while the observed radiation is 
clearly non-thermal, one can nevertheless claim it to be 
related to the optical horizons (blocking regions) via a 
modification of Eq. ([I]). However, this blocking con- 
dition, which at first glance seems rather straightfor- 
ward, in fact contains a number of subtleties once op- 
tical dispersion is introduced. Indeed in situations of 
dispersive propagation with hq = rio(A), where A is now 
the wavelength of the specific mode under considera- 
tion, the simple refractive index blocking condition of 
Eq. becomes a wavelength-dependent "window con- 
dition" defining those modes for which a blocking region 
exists: 
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The window condition ^ will be satisfied by a certain 
set of wavelengths, a set which depends on which form 
for the dispersion relation one chooses as well as on the 
experimental parameters such as the character of the RIP 
itself. However, not all dispersion relations will necessar- 
ily lead to windows — the "window region" might (and 
often does) define an empty set. 

The relevant wavelength dependence of the refractive 
index for fused silica is adequately approximated by the 
Sellmeier relation 
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(Except in the immediate vicinity of resonances where 
the unphysical poles are regulated in real physical sys- 
tems by other effects such as absorption.) In Eq. Q the 



Bi and A^ are a set of coefficients determined by matching 
experimental data for the dispersion. When truncating 
the sum to three poles, the data is best fit with the values 

Ai = 64.25nm, A 2 = 114.00nm, A 3 = 9938.24nm. 

(5) 

which fix the location of the poles in the optical disper- 
sion. 

Typical values of the induced perturbation in the re- 
fractive index rj are of the order of 10~ 3 , and higher, 
increasing with the intensity of the RIP. Calculating the 
vrip relevant for the experiment requires not only the 
knowledge of the refractive index at the typical RIP fre- 
quency, but also the characteristic Bessel beam factor, 
which was measured to be 9b = 6.8° [HI The factor- 
enters the formula for the velocity as a function of refrac- 
tive index as v s i g nai = c/(gn s ig na i) where q = cos(6*s). 



III. PHYSICAL INTERPRETATION 

In Refs. [H E] the window condition of Eq. (|3| was ap- 
plied with the understanding that i>rip was the group 
velocity of the pulse, whereas the the no (A) was taken 
to be the refractive index for phase velocities of modes 
in the fused silica. With this interpretation the window 
was found to be approximately (870, 940) nm. We stress 
that this result is very strongly sensitive to the precise 
numbers used for the relevant quantities in the window 
condition of Eq. ^ due to the relative flatness of the 
function n (X) over the range of wavelengths of interest 
here. (The Cauchy approximation provides insufficient 
accuracy for determining the window condition and it 
is essential to use the Sellmeier approximation for this 
particular purpose.) Such a window appears to be in 
good agreement with the observed spectra, hence sup- 
porting the authors' claim of the Hawking-like nature of 
the observed radiation. Furthermore, the peak shift and 
spectral broadening with intensity arc also claimed to be 
expected in this interpretative framework. In particu- 
lar, the emission bandwidth is predicted to depend on rj 
which in turn is a linear function of the pulse intensity. 

While these aspects of the evidence are surely support- 
ive of a Hawking-like interpretation for the origin of the 
observed radiation, it was immediately noticed by sev- 
eral researchers in the analogue gravity community that 
certain other features of the experiment are problematic 
from this point of view (see e.g. the Comment |10j . and 
Reply [II]). 

We do not wish here to further enter into that on-going 
debate, but we shall instead limit ourselves to stressing 
two points of main concern from our point of view: 

• First of all we find odd the observation that pho- 
tons are emitted at 90 degrees with respect to the 
ostensible horizons set by the window condition |l]) . 
In fact, the window condition only holds for modes 
aligned with the RIP direction of propagation, so 
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one would not expect at all in this model any pho- 
ton production at 90 degrees to the beam axis. 

• Furthermore, we think that some concern has to be 
raised with the application of the window condition 
itself, as it rather oddly mixes the group velocity 
of the RIP with the phase velocities of the vacuum 
modes in a dispersive medium. On the contrary, 
we feel that there are good physics arguments for 
expecting the window condition to be determined 
by only considering group velocities. 

Specifically, let us consider an observer in the rest 
frame of the pulse as it moves down through the glass 
at velocity vrip < c. Then a probe signal chasing the 
RIP from behind will be seen to be approaching with a 
velocity given by the relativistic formula for the combi- 
nation of velocities 
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(6) 



Including the effect of a Bessel beam factor, we have 
vrip = c/(qn mP ). Then 
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From this we see that, in the optical band, zero approach 
velocity (in the RIP frame) is achieved at the location x* 
such that 
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This defines a wavelength-dependent blocking region 
since the probe sees the presence of the RIP as a lo- 
cal increase in the ambient refractive index. Since q < 1, 
and these are interpreted as group refractive indices, only 
faster wave packets (ones with a smaller group refrac- 
tive index, including the Bessel beam factor q) possess 
a horizon. Since n pro b e (A, x) varies between the back- 
ground value no (A) and the value at the peak of the RIP 
no (A) + n, the above condition Eq. ^ is equivalent to 
Eq. j3j. 

However, from a theoretical perspective, use of the rel- 
ativistic combination of velocities formula seems to us 
to be well-motivated only if i>pjp and fprobe are cither 
both group velocities or both phase velocities. Certainly 
vrip is a group velocity since the RIP is a soliton wave 
packet. This rather strongly suggests that in determin- 
ing the window region i> pro be should also be interpreted 
as a group velocity. 

If we now adopt a group velocity interpretation of the 
window condition, we are obliged to consider the group 
velocity refractive index which is defined as 
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where a; is a frequency. Using this formula the window 
condition becomes 
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The problem is that using this (more physically justified) 
formula, and the previously given values for the charac- 
teristic quantities, there is no range of wavelengths for 
which the window conditions is satisfied. As a function 
of 6b the window is non-empty only for 6 b J$ 3° when 
n = 10 -3 and Arip = 1050 nm. 

In summary, we feel that these two facts alone are quite 
sufficient for motivating serious consideration of alterna- 
tive physical explanations for the interesting observations 
reported in Ref . [5J . Specifically, we shall argue here that 
the dynamical Casimir effect (particle production from 
the vacuum due to a time varying external field) is a 
competitive explanation that certainly deserves further 
theoretical and experimental investigation. 

In particular, we conjecture that the rapidly varying 
refractive index in the fused silica excites vacuum modes 
leading to substantial photon creation. This would be 
naturally isotropic in the subspace orthogonal to the di- 
rection of propagation of the RIP, and would not be di- 
rectly dependent on the above mentioned window con- 
dition. In the following discussion we shall present an 
analytically solvable 1+1 model which, in spite of its sim- 
plicity, is able to reproduce many of the salient features 
of the experimentally observed photon emission. 



IV. SETTING UP THE CALCULATION 

As is standard, we simplify the problem to the scalar 
electric field component of the full electromagnetic field, 
with the understanding that the results derived will ap- 
ply to both photon polarizations independently, and will 
contribute half of the total flux expected from a full calcu- 
lation. The equation of motion for the 1 + 1 dimensional 
electric field <b with a time variable refractive index n is 
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In momentum space this reads 

4> k ,rr + c 2 k 2 n 2 (T)(j) k = 0, (12) 
where we have introduced the "conformal time" variable 



T{t) 
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For the refractive index we choose a sech (•) time- 
dependent perturbation on top of momentum-dependent 
background value 

n 2 (t) = n (k) 2 + 2n n (fc) sech 2 (t/i ) + C(n 2 ), (14) 

where n is a small parameter and no(fc) would be a model 
for the refractive index such as the Sellmeier approxima- 
tion. This choice is motivated not only by exact solubility 
of the model, but also for its close approximation to the 
experimental Gaussian profile reported in [S]. 
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Then the equation of motion becomes 

(f>k,rr + c 2 k 2 [nl(k) + 2r) n (k) sech 2 (i(r)/^)] <j> k = 0, 

(15) 

which can immediately be recognized as describing a 
time-independent Schroedinger scattering problem 

2m 



t** + -J? l E ~ V(x)] V = 0, 



with 



V(x) = Vq sech 2 (a;/xo), m 



c 2 k 2 h 2 , 



while 



E = n (k) 2 , 



and 



V = -2 V n (k). 



(16) 



(17) 
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This analogy allows us to directly write down the so- 
lutions to our equation of motion by appealing to well- 
known textbook results. 

The Bogoliubov coefficient for our mode equation 



(151 is related to the transmission coefficient T for the 
scattering problem (161 by \fik\ 2 = 1/T — 1. From the 



literature (see for example Ref. [23] and many references 
therein) we find the transmission probability 
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(20) 



with r^. = ckr n Q (k). The function rises from zero at 
k = to a maximum /3 max at fc max before again decaying 
to zero for large k. Physically this structure reflects the 
transition from low k modes which experience the RIP 
as a sudden, approximately instantaneous influence, to 
high k modes for which the RIP represents an adiabatic 
process. There are also features at the poles (resonances) 
of the refractive index, but here we focus on the region 
between the resonances, which is the region of physical 
interest due to it being the only region where observations 
in the RIP experiments are preformed. 

Furthermore, as we will show below, all of the features 
of the function deriving from the time dependence 
of n(t) occur for values of k for which no(k) is essen- 
tially constant, whence the effects of dispersion can be 
ignored for the purposes of this calculation. Accordingly, 
we henceforth set no(k) — uq ~ 1.458 for what follows. 

Due to the smallness of w we can expand the cosine 
about 7r/2 giving 
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( 21 ) 

where we see explicitly the quadratic dependence of the 
spectrum on r\ at lowest order. 



SPECTRUM AND FLUX 



The 3 + 1 dimensional photon flux density is given by 
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(22) 



Special care should be taken when converting to wave- 
lengths (which is relevant for the comparison with the 
literature) . The peak wavelength of emission A pea k is not 
simply given by A(/c poa k) since A pea k is the maxima of the 
function which enters the integral over A, including the 
Jacobian factor 



dN 



dVo\ d\ 



dk 
dX 



dN 



2ir dN 



dVo\ dk A 2 dVol dk 



k=2ir/\. 



From the lowest order term in Eq. (21) we get 
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Noting that the function l/(u 8 sinh 2 (l/tt)) has its maxi- 
mum at a point extemely close to u = 1/4 we obtain 



x 1 2 f 
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(26) 



Above we have simply bootstrapped our initially 1+1 
dimensional results to 3+1 dimensions using the inter- 
pretation of the equation of motion as applying to the 
magnitude k of the wave vector k. This is common pro- 
cedure in dynamical Casimir effect calculations due to 
the isotropic nature of the emitted radiation. We have 
in mind the idea that the detector which observes the 
emitted photons actually watches only a small element 
of the silica glass which, as the pulse passes through, ap- 
pears to the detector to undergo a rapid time variation in 
its refractive index. Within our simple model this boot- 
strapping can be effectively achieved by integrating our 
Bogoliubov coefficient \Pk\ 2 over a three dimensional k 
space, as we have done above. 

Calculating the total flux we have 
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The emitting region in this 3+1 dimensional interpreta- 
tion would approximately be the cylindrical region de- 
fined by the width of the laser beam and length of the 
region over which the refractive index is rapidly vary- 
ing. Observationally the radius of the laser beam is 
W ~ 10~ 5 m 22 , and we can estimate the length of 
the active region as ctono, so the physical volume of the 
active region is V = Trx^cto "■())■ The total number of 
emitted photons is estimated to be 
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VI. DATA FITTING 

The parameters given explicitly in Ref. [3] are 

n = 1.458, t = 2.5 x 1(T 14 s, r) = 10~ 3 , (29) 

where to is calculated from the reported initial spatial 
size of the pulse at its production xo ~ 10~ 5 m and 
the RIP velocity. Note that, up to a small perturbation, 
r(t) = t/riQ which we use inside the sech 2 (-) function 
to maintain the integrability of the equation of motion. 
This effectively changes the above initial value of to to 
to/nl = 1.2 x 1CT 14 s. 

In principle these factors are all we need to derive our 
estimate of the peak frequency and flux. However, it is 
well known that the RIP shape is not invariant along 
its propagation (see e.g. Ref. [I]). Indeed, there are two 
additional effects which are certainly necessary to take 
into account when fitting the data. These are "pulse 
steepening" and "pulse cresting" which we shall discuss 
below. 



A. Pulse steepening 

Pulse steepening is a non-linear optical effect, akin to 
wave shoaling in water waves, which can be responsible 
for altering the size of the spatial region over which the 
pulse goes from its background value to its peak value 
and back by a factor of ten or more. This effect works 
oppositely on the trailing and leading edges of the RIP; 
while the trailing edge becomes steeper the leading edge 
typically becomes shallower. 

An important question is what, if any, physical mech- 
anism is responsible for setting an ultimate limit to the 
pulse steepening. This point is not completely clear in the 
extant literature, although in many experimental situa- 
tions one adopts the rule of thumb that pulse steepening 
saturates at "about twice the carrier frequency" [H [22] . 

We believe, however, that it is reasonable to expect 
that the limiting mechanism should be more closely re- 
lated to physical properties of the silica glass. Further- 
more we wish here to explore which effects can provide 
absolute limits to the pulse steepening and hence to the 
peak frequency predicted by our model. In this respect 
there are two obvious candidates. One is the plasma fre- 
quency uj p and the other is optical absorption, known 
to be particularly effective near poles in the dispersion 
relation (here characterized in terms of the Sellmeier re- 
lation). 

The plasma frequency definitely provides an ultimate 
limit to pulse steepening, as it is fundamentally the elec- 
trons in the optical glass which communicate the propa- 
gating electromagnetic wave, and which themselves inter- 
act on a time scale given by the plasma frequency. From 
the Sellmeier relation, and the definition of the plasma 
frequency as the coefficient of the leading order deviation 



as the refractive index goes to 1, 

"2 A 2 
n 2 ^l--f = l--, (30) 

U A p 

one obtains 

«* = V^; A * = 7SW (31) 

For the suprasil glass used in the experiment one gets 

Up = 2.6 x 10 16 Hz, A p = 72.7 nm . (32) 

With the above values one gets a typical timescale for the 
steepening of about 2n/u p ~ 2.4 x 10 -16 s. This would 
be the ultimate limit to pulse steepening and we would 
expect the maximum steepness to closely approach, if not 
exactly saturate, this bound. 

The function sech 2 (i/io) varies between its minimum 
and its maximum over a time scale of 2.4 x 10~ 16 s when 
the parameter t is 2 x 10~ 16 s. 

The above estimate represents an absolute upper 
bound to the steepening. However, one can argue that 
before the plasma frequency can start to play any role, 
pulse steepening will be limited by dispersive effects and 
in particular by the pole of the Sellmeier relation at ap- 
proximately 114 nm (see Eq. Q). Such poles represent 
resonances inside the silica glass and their role in limit- 
ing the time variation or propagating signals is essentially 
identical to that of the plasma frequency. Nonetheless, if 
one takes, as we shall do, twice the time scale set by the 
plasma frequency to ~ 4 x 10~ 16 s, the resultant steep- 
ness is just below that set by the Sellmeier pole scale, 
exactly in the right ballpark to model a limiting scale set 
by the pole's presence. Specifically the propagating RIP 
varies between its maximum and minimum over a length 
scale of approximately 140 nm when to ~ 4 x 10~ 16 s 
while the Sellmeier pole sits at 114 nm. 

Let us here anticipate that to match the data our 
model requires that the pulse steepening saturates near 
the above scale to rather than, say, at twice the car- 
rier frequency given set by the above-mentioned rule of 
thumb. This is a relevant point as it might be crucial 
in dismissing or supporting the present proposal against 
others. 

B. Crest amplification 

Another important effect which we expect to occur in 
a propagating and steepening pulse is crest amplification, 
whereby the maximum amplitude of the pulse increases 
as the RIP steepens. The cause of this effect can be 
understood by a conservation of energy argument. First, 
the integral of the intensity is constant (and given by the 
energy of the RIP). Second, due to the non-linear Kerr 
effect the intensity is proportional to the perturbation in 
the refractive index Sn(x). Therefore 

J Sn(x) dx = constant, (33) 



G 



so that 



Casimir spectrum as a function of wavelength 



rj J [shape of pulse] (x) dx — constant. 



(34) 



That is, if the pulse changes shape in such a way as to 
make the area under its normalized curve lower, we can 
expect the maximum amplitude, controlled by the pa- 
rameter rj to increase. Note that the total flux of emitted 
quanta is (approximately) quadratically dependent on rj. 

While in the current context we lack any detailed un- 
derstanding of the size of this effect, we argue that one 
can easily expect 77 to increase up to one order of magni- 
tude. We shall hence use 77 ~ 10~ 2 in what follows. 

In the closely related but distinct context of fibre-optic 
amplifiers, power amplifications of 20 dB, 30 dB, and 40 
dB are not uncommon [Mj-f2|Jj. These are power ampli- 
fications of 10 2 , 10 3 , and 10 4 respectively corresponding 
to amplitude increases of 10, 32, and 100 respectively. So 
a cresting factor of up to 100 (corresponding to rj ~ 10 _1 ) 
is not inconceivable. 

Returning to the current context, the pulse steepening 



-16 



s. 



we have argued for, from to ~ 10~ 14 s to to ~ 10 
will (provided most of the pulse is concentrated in a nar- 
row spike with a long shallow leading edge) squeeze the 
pulse temporally by up to a factor of 100, thereby po- 
tentially increasing the amplitude by a factor of up to 
100. Again, this implies that a cresting factor of up to 
100 (corresponding to 77 ~ 10 _1 ) is not inconceivable. 
Therefore, we argue that our choice of rj ~ 10~ 2 is rather 
conservative. It is important to note that crest amplifi- 
cation mainly affects the total flux, at 0(rj 2 ), and that 
effects on the peak wavelength and width of the spectrum 
are relatively subdominant, at 0(r/ 3 ). 



C. Peak frequency and Flux matching 

We are now in shape for computing the salient features 
predicted by our model. From the previous calculation 
of the peak frequency we get (using the least amount of 
pulse steepening, that is the largest value of f in the 
above given interval) 



1 „ 

Apeak = 2 71 " ct o n o — 858 nm. 



(35) 



This is in striking agreement with the observational data. 

Note again, however, that to get this nice matching for 
the peak frequency it is essential that the pulse steepen- 
ing occur on a timescale and distance scale close to those 
set by the pole in the Sellmeier relation discussed above. 
In fact we have been generous and chosen a time scale 
slightly longer than the Sellmeier resonance time scale. 
This has the side effect of avoiding the resonance at 114 
nm — which is the first resonance one encounters in the 
Sellmeier approximation, (moving up in frequency from 
the carrier frequency as the pulse steepens) . Furthermore 
if one adopts the numerical model discussed in [9], then 
Fig. 10 of that reference implies a timescale of « 1 fs at 
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FIG. 1. Emission spectrum for the dynamical Casimir effect, 
as given by Eq. ( 25 1 and using the parameter values discussed 
in the text: 



n = 1.45, t = 4 x 10~ ib s and 77 = 10 . 
Recall that the total flux is given by the area under this curve 
multiplied by the volume of the emitting region. 



the trailing edge of the RIP; this is within a factor 2| 
of our preferred value. In contrast, if one simply uses 
some small multiple of the carrier frequency as an esti- 
mate of the pulse steepening, the fit is quite bad. It is 
therefore important to develop a clearer understanding 
of the amount of pulse steepening to be expected in this 
particular experimental setup. 

On the other hand the total flux becomes 



N 



8tt 3 rj 2 W 2 
c 2 t\n\ 
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0.6 photons. 



(36) 



This estimate should be supplemented by a factor two to 
take into account the two photon polarizations and with 
an additional solid angle scaling f2/(47r) to account for 
the finite size detector. In the experiment the detector 
was a 5 cm diameter lens placed 10 cm from the fila- 
ment [22 . In terms of the half-angle subtended by the 
detector one has 57 = 27r(l — cos#) whence H/Att w 0.016. 
Unfortunately we do not have an accurate account of the 
observed flux. 

We consider the above numbers to be somewhat en- 
couraging and supportive of the dynamical Casimir effect 
explanation. In Fig.[T]we show the emission spectrum for 
the dynamical Casimir effect as a function of wavelength 
as given by the function F(-) in Eq. (25). 



D. Peak shift and broadening 

Of course, having shown that the relatively unsophis- 
ticated and crude model here proposed can can get any- 
where near reproducing the observed peak frequency, and 
provide anywhere near a reasonable flux, while encour- 
aging, is certainly not enough for claiming to have ex- 
plained the experiment under investigation. In particu- 
lar, we have not yet explained the abovementioncd shift 
of the peak frequency and broadening of the spectrum. 
We shall now argue that our model can account for the 
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Peak wavelength of emission as a function of input energy 
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We note that with an increase of the input energy by 
a factor of 10, our model predicts an increase in flux by a 
factor of (approximately) 100. In the case of a factor of 
5 increase, as in the experimental data, we would expect 
from our model a flux increase of a factor of (approx- 
imately) 25. This predicted behaviour is qualitatively 
consistent with the observed spectra shown in Fig. 3 of 
Ref. [3] where one sees an increase in peak flux from ap- 
proximately 2 photoelectron counts at the lowest input 
energy (270 /iJ) to approximately 50 at the highest input 
energy (1280 fiJ). 



FIG. 2. Shift of the location of the maxima of the emission 
spectrum as a function of wavelength (in nm) for the sech 2 (-) 
profile as the input energy is increased, here modeled by in- 
creasing the parameter n from 2.5 x 1CP 3 to 10 _1 . 



former, and we shall critically discuss the relevance of the 
latter. 

Experimentally, working from Fig. 3 in Ref. [3], we can 
summarize the above two points with the following two 
empirical formulas: 

A pcak ~ 850 nm + 25 nm/mJ x E. (37) 



a x ^ 100 nm/mJ x E. (38) 

where a is the full-width-at-half-maximum (FWHM) of 
the Gaussian fits and E is the input beam energy. 



Peak shift 



From Eq. (21) we see that the peak wavenumber (and 



hence also the peak wavelength of emission) depends on 
rj when we move to higher orders in the 77 expansion of 
the exact result. In Fig. [2] we plot the maxima of the 
spectral function F(X) for a range of values of 77 between 
0.25 x 10 -2 and 10 _1 . Note that this shifting is consistent 
with the shift observed in the Belgiorno et al article [3J; 
increasing the input energy shifts the peak emission wave- 
length to a higher value. We also note that the shifting 
is not exactly linear once we also take into account the 
rf term in the expansion of | (3k \ 2 ■ 

From Fig. 3 of Ref. [3] one might argue that the peak 
wavelengths of emission are approximately 

Apeak ^ 860 nm at I = 270 //J, (39) 
Apeak ^ 875 nm at I = 1280 fiJ. (40) 

That is, the experiments observe a 15 nm shift in peak 
wavelength from a factor of 5 increase of the input energy. 
From our result shown in Fig. [2] we see that our model 
predicts a spectral shift of about 10 nm from a factor of 
10 increase of the input energy (given by 77). 



2. Peak broadening 

Peak broadening, as we have mentioned above, is a 
feature observed in the data by matching the measured 
spectrum to Gaussian curves by eye. It is interesting to 
note that although the magnitude of the effect does not 
seem to be in agreement with the observational data, our 
oversimplified model also predicts a spectral broadening 
as the input energy (77) is increased. We use the FWHM 
definition for bandwidth, and find that a\ behaves as 



-2 



a\ ~ 767 nm for r\ = 10~ 
<t\ ~ 771 nm for 77 = 10 
~ 787 nm for rj = 10 _1 



(41) 



As an increasing function of rj, we see that the predicted 
broadening is in qualitative (though not quantitative) 
agreement with the data. The broadening predicted by 
our simple model is shown in Fig. [3j where we plot the 
scaled spectrum: 



[-P^oMlpeak 



(42) 



That is, we plot the spectrum scaled to the maximum of 
the spectrum at 770 = 10~ 3 , multiplied by an extra factor 
of (770/77) 2 . This is done in order for the spectral curves to 
be visually comparable, and to observe the non-quadratic 
dependence of the spectrum on 77. Note that in this figure 
we have used an exaggerated range of values for the pa- 
rameter 77 to highlight the effect. We also note that this 
figure shows the shifting peak wavelength of emission, 
where we see the peak moving to higher wavelengths as 
77 is increased. We again emphasize that the reported 
broadening in Fig. 3 of Ref. [3] arises from matching 
noisy data to purely phenomenologically chosen Gaus- 
sian curves (there is no physics reason underlying the 
choice of Gaussian curves), and that the reported magni- 
tudes might, in a slightly more conservative analysis, be 
lower than those reported. 



VII. CONCLUSIONS 

The model provided in this article is very simple, in 
many ways overly simple, but nevertheless provides a tol- 



Scaled flux spectra as a function of wavelength 
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FIG. 3. Normalized spectra F V (X) x (vo/v) 2 /{Fvo MLeak- 
Here F(-) is obtained from Eq. ( |25[ ), but modified to include 
the higher order terms indicated in Eq. ( |2I[ ). We normalize 
to the peak of the spectrum for 770 = 10~' 3 , and display the 
result as a function of wavelength in meters. By multiply- 
ing by (?7o/??) 2 , we can see how the shape and maximum of 
the spectrum as a function of 77 is not exactly quadratic in 77, 
but in fact depends on the higher order terms in indicated in 
(21 1. A numerical routine to find the FWHM bandwidths of 
these curves (which are i nde pendent of the amplitude) give 
the results shown in Eq. (41 1. 



erable fit to many of the known aspects of the experimen- 
tal situation. The dynamical Casimir effect envisaged in 
this model has a long and convoluted history in its own 
right (see for instance [2"TH3"3"] V Particularly attractive 
features of the dynamical Casimir model are that it does 
not need any window, that there is no preferred direc- 
tion, and that in performing the experiment one could 
look from any angle. 

In the present context, the most pressing aspect of the 
model is to more fully understand pulse broadening, and 
to undertake a detailed confrontation of the model with 
empirical reality. (Specifically, it would be important to 
develop precise quantitative estimates of the amount of 
pulse steepening and crest amplification to be expected 
in the relevant experimental context.) Another poten- 
tial line of development would be to go to a full 3+1 
dimensional calculation — we have already undertaken a 
perturbative 3+1 dimensional calculation but it taught 
us nothing new beyond the present simple model. 

From a theoretical perspective, one might argue that 
the natural frame of reference for understanding this ex- 
perimental result is the RIP frame. Of course, in this 
frame the simple picture of a time dependent refractive 



index seems to be lacking, but one instead has to view the 
photons as being created by the optical glass rapidly flow- 
ing through the RIP; with the quantum ground state of 
the electromagnetic field being defined by the rest frame 
of the suprasil optical glass, quantum excitations are now 
generated by motion through an inhomogeneous medium. 

Furthermore, in the RIP frame, modes endowed with 
modified dispersion relation will see a sequence of sub- 
luminal regions of different intensity. Situations like this 
have been shown to still lead to particle production even 
in absence of any horizon/blocking region [54T43T] for the 
"dual" situation of modes subject to a supersonic disper- 
sion relation in a "flow" with two supersonic regions of 
different speeds. Indeed there seems to be mathemati- 
cal correspondence that leads to the same behavior when 
passing from super- to sub-sonic in the dispersion rela- 
tion and in the character of the flow regions [38]. If con- 
firmed this analysis would corroborate our claim that a 
particle production from the quantum vacuum could still 
take place in our system even in the absence of analogue 
horizons. 

Finally, we emphasize that the dynamical Casimir ef- 
fect is still quantum vacuum radiation, just not Hawking 
radiation. Both the dynamical Casimir effect and Hawk- 
ing radiation produce squeezed state outputs. Direct 
tests of the difference between these two effects should 
focus on the different correlation structure in the emit- 
ted photons. Specifically, it would be interesting to con- 
sider 2- photon correlations, and back-to-back photons 
(see e.g. [3H I3S]). These effects are likely to be differ- 
ent in detail. One might also need to do considerably 
more than just looking at the squeezed state properties. 
At a very practical level, the output radiation might not 
be "clean" , due for instance to interaction with silica de- 
fects. While much theoretical work remains to be done, in 
many ways the ball is now in the experimentalists court. 
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